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Abstract 

In Kolmogorov's phenomenological theory of turbulence, the energy spectrum in the 
inertial range scales with the wave number k as A; -5 / 3 and extends up to a dissipation 
wave number k u , which is given in terms of the energy dissipation rate e and viscosity v by 
k v oc (e/ z/ 3 ) 1 / 4 . This result leads to Landau's heuristic estimate for the number of degrees of 
freedom that scales as Re 9,/4 , where Re is the Reynolds number. Here we consider the possi- 
bility of establishing a quantitative basis for these results from first principles. In particular, 
we examine the extent to which they can be derived from the three-dimensional Navier- 
Stokes system, making use of Kolmogorov's hypothesis of finite and viscosity-independent 
energy dissipation only. It is found that the Taylor microscale wave number kx (a close 
cousin of k v ) can be expressed in the form kx < CU/u = (CU/ ||it||) 1 / 2 (e/z>' 3 ) 1 / 4 . Here U 
and ||u|| are, respectively, a "microscale" velocity and the root mean square velocity, and 
C < 1 is a dynamical parameter. This result can be seen to be in line with Kolmogorov's 
prediction for k v . Furthermore, it is shown that the minimum number of greatest Lyapunov 
exponents whose sum becomes negative does not exceed Re 9//4 , where Re is defined in terms 
of an average energy dissipation rate, the system length scale, and v. This result is in a 
remarkable agreement with the Landau estimate, up to a presumably slight discrepancy 
between the conventional and the present energy dissipation rates used in the definition of 
Re. 
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1 Introduction 



For the past several decades, Kolmogorov's phenomenological theory of turbulence 1 has 
been a starting point for a majority of theoretical, numerical, and experimental studies of 
fluid motion. A cornerstone of the theory is the notion of energy cascade and dissipation 
anomaly. This means that no matter how small the viscosity, the large-scale kinetic energy 
that generates and sustains the turbulence can be transferred to the correspondingly small 
scales for removal by viscous effects. More precisely, the energy dissipation rate has been 
conjectured to remain finite and nonzero and to become independent of viscosity in the 
limit of small viscosity. This conjecture was first extended by Obukhov 2 and Corrsin 3 to 
the case of passive scalar transport and mixing in turbulent flows, for which the dissipation 
of the scalar variance by molecular diffusion has been thought to be diffusivity indepen- 
dent in the limit of vanishingly small diffusivity. In a similar manner, Batchelor 4 adapted 
Kolmogorov's theory for two-dimensional turbulence, by hypothesizing that the dissipation 
of enstrophy (half mean square vorticity) becomes independent of viscosity (and remains 
nonzero) in the inviscid limit. Like the energy dissipation rate in the Kolmogorov theory, 
the supposedly finite dissipation rates of the scalar variance and enstrophy are key param- 
eters in the Obukhov-Corrsin and Batchelor theories. Since then, the concept of cascade 
and dissipation anomaly has been thought to apply to a variety of fluid systems in other 
contexts. Thus, the Kolmogorov theory of turbulence evidently has become one of the most 
influential theories in science. 

An important prediction of Kolmogorov's theory is that the energy cascading range 
(known as energy inertial range) scales with the wave number k as k~ 5 ^ 3 and extends up to 
a dissipation wave number k u , which is given in terms of the viscosity- independent energy 
dissipation rate e and the viscosity v by k v oc (e/ z/ 3 ) 1 / 4 . This wave number marks the end of 
the energy inertial range, beyond which viscous forces become significant and the dissipa- 
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tion of energy mainly occurs. Heuristically, if the turbulence is homogeneous and isotropic, 
consisting predominantly of vortices of volume k~ 3 , then there are (Lk u ) 3 oc (L 4 / 3 e 1 / 3 /z/) 9 / 4 
such vortices within a domain region of size L. This is the basis for the Landau 5 estimate 
iV oc Re 9//4 for the number of degrees of freedom N, where Re = L 4 / 3 e 1 / 3 /z/ is the Reynolds 
number. A number of authors 6-11 have attempted to address these results from first prin- 
ciples, pending a verdict on solution regularity. The reported upper bounds for N in the 
literature have a wide-ranging dependence on Re, depending on how Re is defined (and to 
some extent on how iV is defined). Constantin et al. 6 (see also Foias et al. 7 ) showed that 
the number of determining modes is proportional to Re 9 / 4 , where e in Re is replaced by 
the asymptotic average of the supremum of the energy dissipation rate. They also found 
that the attractor dimension scales as Re 3 , where Re is defined in terms of the maximum 
fluid velocity. This dimension was found by Gibbon and Titi 9 to scale with the presently 
defined Re as Re 18 / 5 . Gibbon 11 showed that the local number of degrees of freedom scales 
as Re 3 , where Re was defined in terms of the local fluid velocity. Among these results, the 
first one appears to be closest to the Landau estimate, differing from it by the definition of 
e, not by the power of Re. 

In this study, we revisit the Kolmogorov prediction k u oc (e/z/ 3 ) 1 / 4 and Landau esti- 
mate N oc Re 9//4 , by further examining the extent to which these can be deduced from 
the three-dimensional Navier-Stokes system, making use of Kolmogorov's hypothesis of 
finite and viscosity-independent energy dissipation only. It is found that the Taylor mi- 
croscale wave number kr, which can be thought of as a dissipation wave number, satisfies 
kr < CU/v. Here C is a nondimensional dynamical parameter satisfying C < 1 and U 
denotes a "microscale" velocity — a time-dependent velocity scale associated with the en- 
strophy production term. This result is derived without an a priori knowledge of e. Upon 
invoking Kolmogorov's hypothesis, we obtain k T < (CU / ||it||) 1 / 2 (e/i/ 3 ) 1 / 4 , where ||n|| is 
the root mean square velocity, which can be seen to be in line with the prediction for k v . 
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Furthermore, it is shown that along any given bounded trajectory, the minimum number 
of greatest local Lyapunov exponents whose sum becomes negative is bounded from above 
by Re 9//4 , where Re is defined in terms of a new domain- average energy dissipation rate. 
This result means that finite-dimensional volume elements in the infinite-dimensional phase 
space (solution space) contract exponentially if their dimensions exceed Re 9 ^ 4 . This is also 
an upper bound for generalized dimensions, such as the box-counting or Hausdorff dimen- 
sions, of a nontrivial attractor (in forced turbulence) if one exists. This result is consistent 
with the Landau estimate, up to a difference between the usual and present average energy 
dissipation rates used in the definition of Re. This difference, which is arguably slight, is 
the "extent" to which the present analysis needs to reach for a complete agreement with 
the Landau estimate. Thus, the present result is a step closer to that predicted by the 
classical theory. 

2 Preliminaries 

The motion of an incompressible fluid is governed by the Navier-Stokes equations, 

Ut + (u- V)it + Vp = uAu, (1) 
V-u = 0, 

where u(x,t) is the fluid velocity, p is the pressure, and v is the viscosity. We consider 
unforced dynamics for convenience. The results can be seen to carry over to the forced 
case without change, except for a minor modification to section 3, where the forcing term 
would be ignored if the injected enstrophy is negligible compared with that produced by 
the vortex stretching mechanism. This is the case for forcing that injects energy at large 
scales. 12 Equation (1) is considered in a periodic domain D = [0, 2ttL] 3 , and all fields are 
assumed to have zero spatial average. This allows each component of u to be expressible as 
a Fourier series in terms of exp{ik-x/L}. Here k = (k±, k 2 , k 3 ), where k±, k 2 , k 3 are integers 
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not simultaneously zero. The incompressible velocity field u(x, t) can then be represented 
by 

u(x,t) = ^w(M)exp j^-^j , (2) 

where k ■ u(k,t) = for incompressibility and u(k,t) = u*(—k,t) for reality. In other 
words, the solution vector-valued function space (phase space) can be spanned by the 
infinite basis {exp[ik ■ x/L]}j~- Throughout this study we assume strong solutions up to 
the time under consideration. The possibility of subsequent development of singularities 
is not an issue. This allows us to make use of usual assumptions within the realm of the 
classical theory, such as bounded velocity and vorticity. 

The solution space is equipped with the scalar product (it • v) and energy norm ||it|| 
given, respectively, by 

(u-v) = (2nL)- 3 Q^tt- vdxj (3) 

and 

||u|| = (M 2 > 1/2 . (4) 

The advection term (it • V)u conserves the kinetic energy ||n|| 2 /2. More generally, one has 
by integration by parts the identity 

(v ■ (u ■ V)w) = -(w ■ (u • V)w), (5) 

for all admissible v and w. Also by integration by parts followed by the Cauchy-Schwarz 
inequality, we have 

||Vit|| 2 = -(it • Ait) < ||it|| ||Att|| . (6) 

Consider a set of mutually orthonormal functions {v±, v 2 , ■ ■ ■ , v n }, i.e., (vi ■ Vj) = 5ij 
with 8ij being the Kronecker delta symbol. For large n, the number of Fourier modes within 



the wave number radius n 1 ^ 3 / L is approximately n. Their (repeated) eigenvalues under —A 
are k 2 /L 2 = (kf + k% + k 2 )/L 2 < n 2 ^ 3 /L 2 and sum up to approximately n 5//3 /L 2 . Since 
these constitute the smallest eigenvalues of —A, the Rayleigh-Ritz principle implies that 

n „5/3 

El|V^|| 2 >c^, (7) 
where c is a nondimensional constant independent of the orthonormal set {vi, v 2 , ■ ■ ■ , v n }. 

3 Taylor microscale and number of active modes 

In this section we consider the equation governing the evolution of enstrophy || Vu\\ 2 /2, from 
which the enstrophy production and dissipation terms are compared. From this comparison 
during the growing phase of the enstrophy, up to a maximum, we derive two expressions 
for the Taylor microscale wave number and compare them with k u . We then rephrase the 
classical argument leading to the Landau estimate for the number of degrees of freedom. 

By taking the scalar product of the momentum equation of Eq. (1) with Au and inte- 
grating by parts the resulting time derivative term, we obtain 

|| Vu|| y || Vm|| = (Au ■ (u ■ V)w) - v \\Auf . (8) 

The triple-product term in Eq. (8) is responsible for enstrophy production, which holds the 
key to our understanding of turbulence. As it stands, this term is strongly dominated by 
the small scales because of the factors Vu and Au. Indeed, one can readily appreciate this 
claim by rewriting the enstrophy production term in the more familiar form (u> • (u> • V)w), 
where u) = V x u is the vorticity, which clearly involves small-scale quantities only. This 
means that the "macroscale" velocity field plays a largely insignificant role in the enstrophy 
dynamics. In other words, the interactions among the small scales are responsible for intense 
enstrophy production rather than the advection of small-scale eddies and filaments by the 
large-scale flow. This can be true even if the turbulence has energetic large-scale structures 
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superimposed on a sea of small-scale vortices. Motivated by this observation, we define a 
"microscale" velocity U by 

v _ \((Au.(u-V)u)\ _ |(a,.(a,.VH| 
||V«|| ||Aw|| ||V«|| ||A«|| 

As it stands, U can be considered as the typical velocity at the dissipation scale. Since U 
is essentially the domain average of the velocity weighted by two functions of unity norm, 
it satisfies U < HttH^, where HuH^ is the maximum fluid velocity. Given the steep energy 
spectrum of the classical theory, it is likely that U <C H^H^. It may probably be the case 
that U <C \\u\\. Anyhow, with definition (9), Eq. (8) can be rewritten as 

II Vu|| 1| Vu|| < U || Au|| |l V-u.ll - v \\Au\\ 2 



\Auf fUWVuf 



— i/ 1| v«n 2 J . (10) 



" ||Vw|| 2 V \\ Au \ 

During the period of enstrophy increase, up to a local maximum or global maximum, i.e., 
maximum energy dissipation, Eq. (10) implies 

||Vu|| U]Vuf_ CU 
||w|| v ||tt|| ||Aw|| v 

where kx is the Taylor microscale wave number mentioned earlier and C = \\ Vu\\ 2 /(||w|| ||Aw| 
is a nondimensional parameter. By Eq. (6), C satisfies C < 1, but this bound can be highly 
excessive. For example, for the classical spectrum, C tends to zero quite rapidly (as k^ 1 ^) 
in the limit k u — > oo. Note that without a prior knowledge of k v , one could define kx as 
the energy dissipation wave number. Within the context of Kolmogorov's theory, we have 
kr < k v \ and the equality would occur for a k~ x energy spectrum. 

Equation (11) has been derived independently of the Kolmogorov hypothesis of finite 
energy dissipation e. This hypothesis immediately implies that k T = (e/V) 1 / 2 / ||w||. Here, 
we are interested in its dynamical consequences rather than this immediate implication. 
Together with the predicted k~ 5 ^ 3 spectrum, the hypothesis has a profound implication on 



U, i.e., on the enstrophy production term. One can see from Eq. (10) that a finite and 



U <C ||w|| is fully justified within the framework of the classical theory. In any case, Eq. 

(11) can be rewritten in terms of e as 



This brings kx closer to k v in form. Note that if we invoke the Cauchy-Schwarz inequality, 
i.e., replacing C by its upper bound of unity and use U < IM!^, we would have kx < 
(ll^lloo / IMI) 1 ^ 6 /^ 3 ) 1 ^ 4 ' which is apparently an excessive estimate. 

The dissipation wave number k u naturally defines the number of degrees of freedom, 
being the number of dynamically active Fourier modes having wave numbers not exceeding 
k u . This number, denoted by N c , is approximately given by 



where fc = 1/L is the lowest wave number. Equation (13) essentially rephrases Landau's 
estimate mentioned earlier based on the number of vortices having volume k~ 3 within the 
domain D. In what follows, we show that this result largely agrees with an estimate for 
the number of degrees of freedom defined within the context of dynamical systems theory. 

4 Number of degrees of freedom 

This section revisits the notion of number of degrees of freedom from the perspective of 
dynamical systems theory. This number is then estimated for three-dimensional Navier- 
Stokes turbulence, using what is essentially equivalent to the "trace formula," which was 
derived in the 1980s. 13,14 Since then, this formula has become the tool for virtually every 
study on attractor dimension of turbulence in both two and three dimensions. The present 



viscosity-independent energy dissipation rate requires that U || Vtt|| 3 / || Au|| remain finite. 
This implies the scaling U oc k~ l l 3 since |]Vu|| 3 / ||Ait|| oc k\J 3 . Thus, our suggestion that 





(12) 




(13) 
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formulation is a simple version (with other advantages in addition to simplicity, see below) 
of the highly technical formulation that leads to the trace formula. It is an extension of a 
recent study by Tran and Blackbourn 15 on two-dimensional turbulence and is summarized 
in what follows. 

Chaotic dynamics are characterized by the sensitive dependence of solutions on initial 
conditions, resulting in a rapid separation of nearby trajectories — solution "curves" — in 
phase space. In a neighborhood of a given point on a given trajectory (a given solution 
at a given time), this separation is greatest in the orthogonal "directions" corresponding 
to the greatest local Lyapunov exponents. These constitute the most unstable directions 
of the dynamics linearized about the solution under consideration. In general, these di- 
rections and the associated exponents can change continuously along the trajectory, an 
underpinning feature of dynamical complexity. In an infinite-dimensional dissipative sys- 
tem, the number of positive local Lyapunov exponents along a given bounded trajectory 
is presumably finite, followed by a spectrum of negative exponents corresponding to sta- 
ble orthogonal directions. The smallest number of greatest exponents whose sum becomes 
negative (hereafter denoted by N) is significant as phase space n-dimensional volume el- 
ements along the trajectory contract exponentially for n > N. When N is common for 
all points on an arbitrary bounded trajectory, volume contraction becomes universal on 
bounded sets of phase space. This number is an upper bound for the so-called Lyapunov 
or Kaplan and Yorke dimension 16 ' 17 of an attractor if one exists. It is also an upper bound 
for other generalized dimensions, such as the box-counting and Hausdorff dimensions, of 
the attractor. Such an N represents the number of degrees of freedom of the dynamical 
system in question, in the sense that its chaotic dynamics can be adequately described by 
an ^-dimensional model. This makes sense even for cases in which no nontrivial attractors 
are known to exist. Furthermore, N is well defined regardless of whether or not conven- 
tional Lyapunov exponents exist. Another advantage of the present formulation is that the 
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problem of global regularity of solutions of the Navier-Stokes system is not an issue. The 
reason is that N is determined pointwise in time, therefore remaining valid up to the time 
of solution blowup should this turn out to be the case. 

4.1 Local Lyapunov exponents 

Consider the linear evolution of a disturbance v to the solution u(x,t) of Eq. (1) com- 
mencing from a smooth initial field u = u(x, 0). The governing equations for v are 

v t + (u ■ V)w + • V)n, + Vp' = vAv, (14) 

V • u = = V • v, 

where p' is the perturbed pressure. By taking the scalar product of v with Eq. (14) and 
noting that both (v • (u • V)u) and (v ■ Vp') vanish, we obtain the equation governing the 
evolution of \\v\\, 

\\ v \\ ~t~ \\ v \\ — —{v ■ (v ■ V)u) — v ||Vv|| 2 . (15) 
at 

It follows that 

A = ^- In || v || = Uv ■ (v ■ V)«) + v ||Vu|| 2 ) . (16) 
at \\v\\ ^ ' 

Here A is the exponential growth or decay rate of ||u||. 

The greatest local Lyapunov exponent and the corresponding most unstable direction 
can be found by maximizing A with respect to all admissible disturbances v. We denote by 
(Ai, Vi) the solution of this problem, where for convenience (and without loss of generality) 
Vi has been normalized, i.e., = 1. The second greatest exponent A2 and the corre- 
sponding second most unstable direction v 2 orthogonal to V\ can be obtained by maximizing 
A with respect to all disturbances v subject to the orthogonality constraint (v • Vi) = 0. 
Similarly, the pair of third greatest exponent A3 and third most unstable direction v% can 
be obtained by solving the same maximization problem, where the admissible disturbances 
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v satisfy the constraint (v.v±) = (v.v 2 ) = 0. By repeating this procedure n times, we 
obtain the set {vi,v 2 , ■ ■ ■ ,v n } of mutually orthonormal functions and the corresponding 
ordered set of exponents Ai > A 2 > • • • > X n . These can be described more formally by 



Xj = max{ — (v ■ (v ■ V)u) — v ||V«|| 2 } 



V =1 



(17) 



for 1 < j < n, where the maximization is subject to the constraint (v ■ Vi) =0 for 
i < j. These exponents provide a complete picture of solution stability with respect to 
disturbances. In passing, it is worth mentioning that in the maximization problem, the 
solutions Vj arise as compromises between the triple-product and viscous dissipation terms. 
That means that they do not necessarily maximize the former. It would be interesting to 
have a knowledge of the maximizers (A^, v'j) of the triple-product term alone. A comparison 



between the two Lyapunov spectra Xj and A^- and between ||Vvj|| and 
by numerical methods, could provide some invaluable dynamical insights 



, conceivably 



4.2 Upper bounds for the number of degrees of freedom 

We now calculate the number of degrees of freedom N described earlier by minimizing n 
such that the sum Yl]=i Aj is negative. By Eq. (17), we have 

n n 

E A,- = - E (<«; • («; • V)u> + u\\W 3 \\ 2 ) 



E((«-(^-vk>-^||v^|| 2 ), 



(18) 



where Eq. (5) has been used. Similar to the definition of U in the preceding section, we 

define two average quantities U' and Q, respectively, by 

1 



U' 



(n£?=il|V^ 



1/2 



E( M ' ( v i ■ V ) v i) 



(19) 



and 



n = i 

n 



(20) 
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Like U in Eq. (9), U' < juW^ represents a small-scale velocity by virtue of its very definition. 
The reason is that the orthonormal set {vi, i>2, • • • , v n } consists of progressively smaller- 
scale functions Vj, i.e., increasingly greater ||VVj|| as the index j increases. Note that 
a suitable rearrangement of the set may be necessary if it is not already in that order. 
Furthermore, similar to the j-th eigenvalue of —A, we have ||V«j|| 2 oc j 2 / 3 . This implies 
that U' is more strongly weighted by the smaller-scale v/s. It is, however, not known with 
precision how U' and U compare. Now for definition (20), £1 can be thought of as the domain 
average of |Vw| weighted by YJj=\ \ v j\ 2 / n - By definition, the inequality fl < HVitH^ holds. 
On physical grounds, this bound can become excessive for large Re because intense velocity 
gradients are known to be highly concentrated in space, 18 ' 19 effectively getting "moderated" 
under the spatial average in the definition of Q. Moreover, unless the spatial distribution of 
a majority of v/s is strongly correlated to that of Vu (i.e., locally peaking in the same small 
regions as Vu), this moderation can be more effective than that in ||Vw|| = {\S7u\ 2 ) l l 2 . 
The reason is that \Vu\ 2 is more resistant to such moderation than |Vw| as is reflected 
in the fact that (]Vw|) < ||Vw||. Hence, even though Q is undetermined, it is expected 
to be closer to ||Vw|| rather than to HVitH^. This "conjecture" could be readily tested 
numerically, given the linear and kinematic nature of the maximization problem. On an 
optimistic note, it is worth mentioning that one cannot rule out the possibility Q < ||Vw||, 
even though that might seem unlikely. 

Upon substituting Eqs. (19) and (7) into the second equation of Eq. (18) we obtain 




|2 
3 ll 



L • (21 > 
The condition Y^=i ^ requires a straightforward lower bound for n, from which we 
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deduce the bound 

N < c" 3/2 ( — ) = Re 3 , (22) 



'U'L X 3 



where c~ 3//2 has been incorporated into the newly defined Reynolds number. Similar results 
have been reported by Constantin et al. 6 and Gibbon, 11 where their Reynolds numbers 
were defined in terms of H^H^ and of the local velocity \u(x,t)\, respectively. The Gibbon 
estimate (local number of degrees of freedom) becomes the Constantin estimate where 
\u(x, t)\ peaks. 

The above estimate has made no use of the assumption of finite energy dissipation. Now 
if we identify uVL 2 with e, then upon substituting Eqs. (20) and (7) into the first equation 
of Eq. (18), we obtain 

n n 

Xj < nQ — v^2i ll^ill 2 

3=1 3=1 

/ e l/2 cn 2/3X 

^{^-"nry (23) 

The condition X)™ =1 < requires a straightforward lower bound for n, from which we 
deduce the bound 

/ j 4/3,1/3 \ 9 / 4 

JV< C ^ =Re 9/4 ? (24) 

where again c~ 3 / 2 has been incorporated into the Reynolds number Re. This result differs 
from the Landau estimate by the use of fl instead of ||Vw|| in the definition of Re. This 
difference can be slight as argued above. In the classical picture of homogeneous turbulence, 
there would hardly be any distinction between ||Vit|| and Q. 

In passing, it is worth mentioning that the sum of the triple-product terms in the 
first equation of Eq. (18) is quite susceptible to sophisticated (and potentially excessive) 
estimates, which we have thus far deliberately avoided. Consider, for example, the Lieb- 
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Thirring inequality 20 ' 21 concerning the orthonormal set {v\, v 2 , • • • , v n }, 

(n \ 3 / 4 

<c'L 3 / 2 EHV^-II 2 , (25) 



where c' is a nondimensional constant independent of the set {vi, v 2 , ■ ■ ■ , v n }. By applying 
Eq. (25) to the first equation of Eq. (18), via the intermediate step 



5>,- • v)«) < 

J'=l 



|V«|| , 



(26) 



we would arrive at 



£A,<c'L 3 / 2 £|| Vt ,.| 

( 3/4 

<(Ellv«i 



3/4 

||Vu|| — ^E ll^^ill 
j'=i 

' , r3/2 e 1 / 2 c^n 5 ^ 



v 



1/2 



L 1 / 2 



(27) 



where Eq. (7) has been used in the second step. It follows that 



f r li\ 3/5 /r 4/3,1/3 \ 18 /5 



(28) 



where the constant prefactor has been absorbed into Re, which is now defined in terms of 
||Vw|| instead of f2. The price for this is the scaling Re 18//5 instead of Re 9 ^ 4 for N. This 
result (given in a quite different form) was derived by Gibbon and Titi 9 as an upper bound 
for the attractor dimension. 



4.3 Discussion 



In two-dimensional turbulence, N has been found to satisfy 15 



N < C'Re(l + lnRe) 1/3 , 



(29) 



where C is an absolute constant and the Reynolds number Re is defined in terms of the 
materially conserved vorticity, the domain size, and v. Apart from the difference in the 
level of rigor in the definition of Re, there is a sharp contrast between the nearly linear 
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scaling of N with Re in two-dimensional turbulence and the highly superlinear scaling of 
N with Re in the present case. This is due to fundamental differences between the two 
cases. We discuss two most apparent discrepancies in what follows. 

One of these is due to the dimension of the physical space and is easy to recognize. 
Given n Fourier modes of lowest wave numbers, the sum of their eigenvalues under — A— 
a collective measure of viscous dissipation strength — are oc n 5//3 and oc n 2 , in three and 
two dimensions, respectively. This means that for the same Reynolds number, three- 
dimensional turbulence is expected to have more dynamically active modes than its two- 
dimensional counterpart. This makes an intuitively obvious contribution to the difference 
between Eqs. (24) and (29). 

The other contributing factor can be attributed to the discrepancy in the "effective 
degree" of nonlinearity of the small-scale dynamics of the Navier-Stokes equations in these 
cases. In three dimensions, the dynamics are highly nonlinear, effectively quadratic. In prin- 
ciple, the vortex stretching term (w- V)u can give rise to an explosive vorticity growth. 22,23 
On the contrary, the two-dimensional Navier-Stokes system is effectively nearly linear, 
rendering far less intense dynamics of the small scales — a widely recognized fact. One can 
readily appreciate this claim by a quick inspection of the equation governing the vorticity 
gradient Vu;, 

Vu t + (u ■ V) Vcj = uon x Vuo - {Vuj ■ V)u + i/AVw, (30) 
V- « = 0, 

where n is the normal to the fluid domain. In Eq. (30), the sole effect of the first term on 
the right-hand side is to rotate Vcj without changing its magnitude, and the second term 
alone is responsible for vorticity gradient amplification. By ignoring the viscous term for 
convenience, we can deduce from Eq. (30) the equation 

\Vu\ t + u ■ VIVcjI = -t5^t • (Vw • V)« < iVitllVwI. (31) 

I Vcu | 
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This means that following the fluid motion, |Vu;| can grow no more rapidly than expo- 
nentially in time at the instantaneous rate |Vtt|. Now, since vorticity is conserved in the 
inviscid dynamics, |Vw| is relatively well behaved because ||Vw|| = is conserved. In- 
deed, numerical evidence shows that for an initial vorticity reservoir at large scales, |Vw| 
remains largely unchanged up to and beyond the instance of peak enstrophy dissipation. 24 
More precisely, the ratio of the irrotational strain to initially at ~ 2, has been found 
to remain within the range [2,3.5] throughout the said period, during which ||Vu;|| grows 
approximately exponentially by several orders of magnitude. Hence, the vorticity gradient 
stretching term can be said to be marginally nonlinear, if it is to be considered nonlinear at 
all. The same remark can be made about a broad family of fluid systems in the geophysical 
context. For this case, the gradient Vg of the materially conserved potential vorticity q is 
governed by Eqs. (30) and (31), with q replacing u. For this family, the velocity gradient 
Vu is also well behaved. In fact, it is presumably better behaved than its counterpart 
in two-dimensional Navier-Stokes turbulence because ||Vu|| < \\q\\. Thus, the small-scale 
dynamics of this family are effectively marginally nonlinear. 

The enstrophy production in three-dimensional turbulence is a fundamental problem in 
fluid mechanics and has always been a centre of attraction for the turbulence community. 25-3 ' 
This is a formidable problem, being virtually intractable as we have come to realize. In 
the limit of large Reynolds number, analytic and dynamically independent upper bounds 
for (Au • (it ■ V)u) tend to become so excessive that they render no practical value. The 
reason behind these excessive estimates is that in order to bound the norm of a quan- 
tity, say Vu, one usually resorts to norms of its derivatives, such as Au. As a result, 
when ||Vw|| — > oo, its upper bounds usually diverge far more rapidly. Known inequalities 
applicable to the enstrophy production term invariably reduce to the form "1 < oo" (or 
equivalently "0 < 1") as ||Vw|| — > oo. An example is the Cauchy-Schwarz inequality (6), 
which was given in section 2 and briefly discussed in section 3. There, we bypassed this 
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inequality by introducing the dynamical parameter C. Another example is the Agmond 35 ' 36 
inequality WuW^ < c"L 3 ^ 2 || VwH 1 ^ 2 ||Au|| 1//2 , where c" is a nondimensional constant. This 
could be one of the most generous estimates in the present context. Given the excessive 
nature of the available inequalities (when applied to high- Reynolds number turbulence), it 
is desirable, if not crucial, to develop new techniques that could derive dynamically binding 
estimates from the governing equations. 

5 Conclusion 

This study has examined a possible route toward a quantitative basis in support of the 
Kolmogorov prediction for the viscous dissipation wave number k v in three-dimensional 
Navier-Stokes turbulence and the associated Landau estimate for its number of degrees of 
freedom N. For k v , we have taken an indirect approach, by estimating the Taylor microscale 
wave number kr, which is a close cousin of k u . It has been found that kr < CU/u, where U 
is a "microscale" velocity and C < 1 is a dynamical parameter. When expressed in terms 
of the energy dissipation rate e, this result becomes kx < (CU/ ||ti||) 1 / 2 (e/i/ 3 ) 1 / 4 , where ||w|| 
is the root mean square velocity. The latter can be seen to be in line with Kolmogorov's 
prediction for k v . For N, we have taken a direct approach, by deriving an upper bound for 
the minimum number of greatest local Lyapunov exponents whose sum becomes negative. 
The calculations of this bound have been carried out in a general manner, therefore the 
obtained result is universal for bounded trajectories. It is an upper bound for generalized 
dimensions, such as the box-counting and Haussdorff dimensions, of a nontrivial attractor 
(for the forced case) if one exists. It has been found that N satisfies N < Re 9 / 4 , where 
Re is defined in terms of an average energy dissipation rate, the system length scale, and 
v. This result is in a remarkable agreement with the Landau estimate if one identifies the 
conventional energy dissipation rate e = i/||Vt*|| with the newly defined rate vVL , where 
fl is effectively the spatial average of the velocity gradient |Vw| weighted by the average 
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of the squares of mutually orthonormal functions. Although Q is essentially undetermined, 
we have argued that it can be close to ||Vw||. This is the "extent" to which the present 
analysis needs to reach for a complete agreement with the Landau estimate of the number 
of degrees of freedom on the basis of the Kolmogorov theory. In the classical picture of 
homogeneous and isotropic turbulence, there would be virtually no distinction between 
||Vw|| and fi. 
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